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Abstract. We present adaptive optics spectro-imaging observations of the RW Aur jet in optical forbidden lines, 
at an angular resolution of 0.4". Comparison with HST data taken 2 years later shows that proper motions 
in the blueshifted and redshifted lobes are in the same ratio as their radial velocities, a direct proof that the 
velocity asymmetry in this jet is real and not an emissivity effect. The inferred jet inclination to the line of sight 
is i = 46 ± 3°. The inner knot spacing appears best explained by time variability with at least two modes: one 
irregular and asymmetric (possibly random) on timescales of <3-10 yr, and another more regular with ~ 20 yr 
period. We also report indirect evidence for correlated velocity and excitation gradients in the redshifted lobe, 
possibly related to the blue/red velocity and brightness asymmetry in this system. 
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1. Introduction 

Jets from young stars are believed to be powered by ac- 
cretion and accelerated and coUimated by magnetic forces 
(see e.g. Konigl & Pudritz 2000 for a review). However, 
a clear identification of the magnetic launch configura- 
tion is still needed to fully understand the role of jets 
in the physics of accreting protoplanetary disks. The opti- 
cally revealed innermost jet regions of T Tauri stars (TTS) 
hold essential clues to this question, and are the subject 
of growing observational efforts at sub-arcsecond resolu- 
tion (e.g. Lavalley-Fouquet et al. 20021 Dougados et al. 
120001 (hereafter Paper I); Bacciotti, Mundt, Ray et al. 
I2000|l . HH 229, the bipolar jet from the actively accret- 
ing star RW Aur A, is a particularly interesting target. 
It was first identified in long-slit spectroscopy by Hirth 
et al. |[T594|i) . who noted an asymmetry in radial velocity 
between the blueshifted and redshifted lobes of a factor 
~1.8. Further study with HST/STIS (Woitas et al. \MM 
hereafter Paper II) traces this asymmetry down to 0.1" 
« 14 AU from the source, suggesting that it might be an 
intrinsic property of the jet driving engine. In addition, 
the RW Aur jet contains many emission knots, traced as 
far as 50"-100" (Mundt & Eisloffel [TUMjl down to only 
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0.1"-10" from the star (Paper I, Paper II), which may be 
used to test models of knot formation in stellar jets. 

In this Letter, we present a ground-based spectro- 
imaging study of the RW Aur jet that brings several new 
results, including improved knot proper motions, jet incli- 
nation, and knot dynamical times, as well as indirect ev- 
idence for internal velocity and excitation gradients. Our 
observations are summarized in Sect.|21and our results are 
discussed in Sect. O 

2. Observations 

Observations were carried out at the CFHT on December 
18-20, 1998 with the integral field spectrograph OASIS. 
We used a spatial sampling of 0.16"/lens and achieved 
an angular resolution of 0.4" FWHM after adaptative op- 
tics correction. Two spectral ranges were observed (6213- 
6534 A and 6496-6824 A) with a resolving power « 90 km 
s~^. A second dataset with 0.11"/lens was obtained over 
a smaller field of view. Following standard OASIS reduc- 
tion (cf. Lavalley-Fouquet et al. I2000|l . a specific proce- 
dure was applied to subtract at each position the strong 
continuum+cmission spectrum of RW Aur A, scaled by 
the PSF, and to construct datacubes of residual forbid- 
den line profiles. More details are given in Dougados et 
al. (2003, in preparation). Comparison between the two 
lens samplings shows that continuum-subtracted emission 
maps at 0.16"/lens are reliable beyond 0.35" from the star, 
while the 0.11"/lens data are more reliable at d < 0.35". 
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Proper motions in the RW Aur redshifted jet 
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Fig. 1. Illustration of knot proper motions in the in- 
ner redshifted jet of RW Aur. From top to bottom: [S 
II] channel maps obtained in Dec. 2000 with STIS/HST 
(Paper II); integrated [S II] image obtained in Dec. 1998 
with OASIS at 0.16"/lens (this paper); deconvolved [O 
I] image obtained in Dec. 1997 with PUEO at CFHT 
(Paper I). 



Radial velocities are computed in the stellar reference 
frame (heliocentric velocity +16 km s~^; Petrov, Gahm, 
Gameiro et al. I2001II . We adopt hereafter the notation 
[OI]ee[OI]A6300, [SII] = [SII]A6731, and [NII] = [NII]A6583. 

3. Results 



3.1. Proper motions of redshifted knots 

Figure ^ compares our 1998 [S II] image of the red- 
shifted jet with images obtained in Dec. 2000 (Paper II) 
and in Dec. 1997 (Paper I)^. A smooth systematic shift is 
seen between the 3 epochs, indicative of proper motions. 

^ In order to replicate the correct separation between 
RW Aur A and B, the Dec. 1997 image published in Paper I had 
to be corrected for a small 5% underestimate in the FOCAM 
pixel size 
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Fig. 2. Comparison of [S II] images of the RW Aur 
biueshifted jet obtained with OASIS in Dec. 1998 (this 
paper) and with STIS/HST in Dec. 2000 (Paper II), illus- 
trating the proper motion. 

Tabled lists the knot positions and the resulting proper 
motions in arcsec per year. As the OASIS spatial reso- 
lution is not sufficient to detect possible offsets between 
different lines, we adopt for the nominal 1998 position 
of each knot in Table 1 an average of independent esti- 
mates using the strongest two lines ([01] and [SII]) and 
the two lens samplings (0.11" and 0.16"). Knots are num- 
bered from the outside in, starting at the most distant 
one detected by Dougados et al. (Paper I). For the outer- 
most knot seen with OASIS (knot R4), the 3 epochs are 
consistent with a constant proper motion of 0.16" yr~^. 
The proper motion for the following knot (R5) is more 
uncertain, as it appears double in the 2000 STIS image. 
The range of possible values is 0.17"-0.27" yr~^. For knot 
R6, we find again 0.16" yr-\ based on the 1998 and 2000 
epochs only. We note that an innermost emission peak 
is detected in the 0.11"/lens dataset at ~ 0.28". It may 
correspond to the peak at 0.25" in the STIS image, and 
would then be stationary. 



3.2. Proper motions of biueshifted knots 

Figure |21 compares [S II] images of the inner biueshifted 
jet in 1998 (this paper) and 2000 (Paper II). We lack a 



Lpez-Martm et al.: Proper motions in the RW Aur jet 



3 



Knot 




dgs 


doo 


IKadl 




i 






(") 


(") 


(") 


(km s"^) 


(" yr-') 


n 


(yr) 


Red lobe 






<0.28 


0.25 


100-56 


stationary? 






R6 




0.51 


0.82 


M-117 


0.16 


44±8 


1.3 


R5 


1.07: 


1.32 


1.72/1.91 


107-133 


0.24±0.05 


53±9 


8/5 


R4 


2.92 


3.12 


3.45 


97-100 


0.16 


47±5 


18 


R3 


3.9-4.6 








[0.16] 




24-29 


R2 


7.9 








[0.16] 




49 


Rl 


11.1 








[0.16] 




69 


Blue lobe 






<0.19 


0.2 


175-168 


stationary? 






B4 






0.57 


179 


[0.26] 




0.2 


B3 




0.73±0.05 


1.25 


178-168 


0.26±0.035 


45±5 


3 


B2 


2.8 








[0.26] 




12 


Bl 


7.1 








[0.26] 




28 



Table 1. Properties of knots in the RW Aur jet. dgy, dgg, and doo are the projected knot offsets in Dec. 1997, 1998, and 
2000 (Paper I; this paper; Paper II), with respective uncertainties ±0.01", ±0.02", and ±0.05", unless when quoted. 
K-ad is the knot [S II] centroid velocity from this paper (roman font) and from Paper II (italics); ^ the proper motion 
(typical uncertainty ±0.027"yr~-^; values between brackets are assumed); i the angle to the line of sight (error bars 
include uncertainties in both ii and Kad); tgs — tej the time between knot "ejection" and Dec. 20, 1998. 



1997 epoch, as deconvolution residuals in the PUEO image 
dominate over the faint blue jet emission out to ~ 2" 
(Paper I). We find a proper motion of 0.26" yr~^ for the 
main knot (B3). A new knot (B4) has appeared at 0.6" 
in the 2000 STIS image. With a similar proper motion 
to B3, it would have been within 0.05" of the star at the 
date of the OASIS observations and thus unresolvable. 
Finally, our 0.11"/lens dataset shows an innermost peak 
at ~ 0.19", which may correspond to the 0.2" peak seen 
in the STIS image and would then be stationary. 

3.3. Revised jet inclination 

Assuming that knot proper motions trace actual fluid mo- 
tions and not wave pattern speeds, the angle of inclina- 
tion i of the jet with respect to the line of sight is given by 
tan(i) = (/x/"yr~^) x (666.7/T4ad), where V^ad is the veloc- 
ity projected on the line of sight, and a distance of 140 pc 
to Taurus-Auriga is assumed. We adopt as V^-ad for each 
knot the centroid of a gaussian fit to the [S II] line profile. 
Values from this paper and from Paper II (in italics) are 
given in Col. 5 of Tabled Discrepancies in V^ad between 
the two datasets are discussed in more detail in Section 
13.51 The resulting values of i are given in Col. 7, with an 
uncertainty interval that includes both the uncertainty in 
Vrad and that in /i (see notes to Table 1). Inclination val- 
ues are remarkably consistent between knots, within the 
estimated errors. In particular, we derive consistent incli- 
nations for the blue and red lobes, i.e. the proper motion 
in the blue jet is higher than in the red jet by the same fac- 
tor (~ 1.6) as their average centroid radial velocities. This 
fact argues strongly for a real velocity difference between 
the blueshifted and redshifted jets, rather than an appar- 
ent difference due solely to slow bowshock wings dominat- 
ing the emission around the redshifted jet (proper motions 



should then be similar in the two lobes despite differing ra- 
dial velocities) . A weighted average over both lobes yields^ 
i = 46°±3°. 

3.4. Knot spacing, dynamical times, and origin 

The RW Aur observations allow to measure knot spacings 
and dynamical timescales very close to the source and thus 
bring interesting constraints on the origin of knots in stel- 
lar jets. To increase the distance and time span in our 
analysis, we include in Table 1 more distant knots imaged 
with PUEO in 1997, which lie beyond the OASIS field of 
view and were not considered by Woitas et al. (Paper II). 

Table H shows that the knot spacing first increases 
with distance from the star, from 0.7"-0.8" at d ~ 1.5" 
(knots R5-R6, B3-B4) to 1.8"-2.4" at d ~ 3" (knots R4- 
R5, B2-B3), and then levels off at Aa; ~ 3" - 4" for 
d > 4". Adopting as the jet characteristic radius rj half 
the FWHM measured at the same distances in Paper II, 
the inner knot spacing is then Aa; ~ 8 — lOrj. For com- 
parison, the axisymmetric pinch instabilities in MHD jets 
simulated by Lery et al. (2000) produce jet knot spacings 
of Ax/rj < 3, while Kelvin- Hclmholtz instabilities predict 
Ax/rj ~ ^/cs > 10. An even more important observa- 
tional constraint is that the first moving knot appears at 
a distance d„iin < Ax, while the growth scale of an in- 
stability is expected to be several times its wavelength. 
Hence, neither kind of instability seems adequate to ex- 
plain the innermost knot spacings, although a range of 
MHD jet configurations should be explored before they 
are fully ruled out. 

^ The larger inclination value i — 53° found by Woitas et al. 
(Paper II) for knot R4 stems from a proper motion estimate 
20% too high, inferred from the published PUEO image of 
Paper I before correction for pixel scale (see footnote 1) 
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Another possible origin for the knots, which can eas- 
ily yield dmin < Ax, is time variability in the jet ejection 
speed (cf. Equ. (3) in Raga et al. 1990). The last column 
of Tablenigives the dynamical ages of each knot assuming 
constant proper motion (with t = arbitrarily set at our 
OASIS observation date of Dec. 19, 1998). For knots lack- 
ing proper motion measurements, we adopt the typical fi 
in the same lobe (in brackets in Col. 6 of Tabled]). We find 
that the timescale between knots initially increases with 
distance, from 3 to 10 yr, and then appears to level off at 
~ 20 yr beyond 4". Furthermore, the inner knots do not 
appear to be "synchronized" between the two jet sides, 
while there is tentative evidence of synchronization when 
the 20 yr period appears (knots R3 and Bl). We conclude 
that there would be two distinct processes of time vari- 
ability in RW Aur: one that is irregular and asymmetric 
— possibly random (random time variability of the jet 
speed produces asymptotically an increase in knot spac- 
ing with distance; R,aga ll992|l and another more regular 
and symmetrical with a 20 yr period — possibly related 
to a stellar magnetic cycle. We suggest that asynchronous 
variability might occur, e.g., through local magnetic re- 
connexion events that affect only one jet side. 

The apparently stationary knots at 0.2"-0.3" ~ 30-45 
AU from the star need to have a different origin, possi- 
bly recollimation shocks (e.g. Gomez de Castro & Pudritz 
I1993|l or temperature gradients in a progressively heated 
wind (e.g. Garcia et al. l2001|l . but predicted distances are 
presently too model-dependent to yield useful constraints. 

3.5. Velocity gradients and velocity asymmetry 

Figure El plots the [S II] centroid radial velocity as a func- 
tion of distance from the star for the OASIS and STIS 
datasets (the latter is shifted to correct for proper mo- 
tions). In the blue jet, the various datasets agree to within 
10 km s""'^, comparable to expected instrumental and cal- 
ibration errors. The situation is different in the redshifted 
jet: the 0.11"/lens and 0.16"/lens datasets still agree with 
each other within < 10 km s~^, but they lie typically 20 
km s~^ below the STIS values out to the position of knot 
R5. Further out, however, at the position of knot R4, all 
data converge to a very similar T^ad- 

These discrepancies can only be understood if the 
faster gas seen on the redshifted jet axis by STIS (0.1" 
wide slit) suffers from beam-smearing in the OASIS data. 
In other words, sharp transverse velocity gradients must 
occur on scales < 0.4" within the redshifted jet, out to 
distances < 3" from the star^. We note that [N II] cen- 
troids from OASIS (shown as crosses in Fig.OJ agree better 
with the [S II] STIS centroids, despite an angular resolu- 
tion identical to the [S II] OASIS data. This fact further 
indicates that the faster axial gas contributes a greater 
fraction of the total flux in [N II] than in [S II], i.e. that it 
has higher excitation. Interestingly, we find no evidence for 
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such gradients seem confirmed by prehminary analysis of 
transverse STIS slits (Woitas et al., work in progress) 
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Fig. 3. Centroid radial velocity in [S II] (from gaus- 
sian fits) as a function of distance from the star. Filled 
and open circles refer to the 0.1 6" /pixel and 0.11"/pixel 
OASIS data, respectively, integrated over a 1" wide 
slit. For comparison, crosses show centroids in [N II] 
(0.16"/lens), and triangles show [S II] centroids in a 0.1" 
slit with HST/STIS (Woitas et al. 2002), shifted to correct 
for proper motion (by -0.32" in the red lobe, -0.52" in the 
blue; cf. Table 1). 

such velocity and excitation gradients in the blueshifted 
jet. If confirmed, this difference might be related to the 
blue/red velocity asymmetry between the two jet lobes; 
e.g., bowshock wings or turbulent boundary-layers could 
be somehow enhanced in the slower lobe (thus also ex- 
plaining its higher brightness) . Further theoretical and ob- 
servational work is needed to explore this possibility. 
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